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Abstract. It has been recently suggested that a totally asymmetric exclusion process 
with two species on an open chain could exhibit spontaneous symmetry breaking in some 
range of the parameters defining its dynamics. The symmetry breaking is manifested 
by the existence of a phase in which the densities of the two species are not equal. In 
order to provide a more rigorous basis to these observations we consider the limit of the 
process when the rate at which particles leave the system goes to zero. In this limit 
the process reduces to a biased random walk in the positive quarter plane, with specific 
boundary conditions. The stationary probability measure of the position of the walker in 
the plane is shown to be concentrated around two symmetrically located points, one on 
each axis, corresponding to the fact that the system is typically in one of the two states of 
broken symmetry in the exclusion process. We compute the average time for the walker 
to traverse the quarter plane from one axis to the other, which corresponds to the average 
time separating two flips between states of broken symmetry in the exclusion process. This 
time is shown to diverge exponentially with the size of the chain. 
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1 Introduction 



It is well known that systems in thermal equilibrium do not exhibit spontaneous sym- 
metry breaking in one dimension at finite temperature, provided that the interactions are 
short range and that the local variable which describes the microscopic state of the sys- 
tem takes only a finite number of possible values [1] . Common examples of such systems 
are Ising or Potts-like models. On the other hand, when the local variable takes one of 
an infinite number of possible values the model may exhibit phase transitions and sym- 
metry breaking in one dimension. Examples are solid-on-solid models which are used to 
describe wetting phenomena or unbinding transitions of a one dimensional interface from 
an attractive wall [2] . 

Systems which are not in thermal equilibrium, but rather evolve under some stochastic 
dynamics, are different. In many cases these systems reach a steady state which does not 
obey detailed balance. The question is whether such a steady state of a one dimensional 
system can break the symmetry of the dynamical rules under which the system evolves. 
A related problem has been considered in the context of error correcting computation 
algorithms, and an example of a one dimensional array of probabilistic cellular automata 
which exhibits nonergodicity (as would a model with spontaneous symmetry breaking) 
has been constructed [3]. However, this example is rather complicated and not widely 
understood. 

Recently a simple model was introduced [4, 5] , describing a totally asymmetric exclusion 
process for two species of particles moving on an open chain of length N. Each lattice 
site may be occupied by either a positive (-I-) or a negative (— ) particle, or by a hole (0). 
Positive particles move to the right, negative particles to the left. During an infinitesimal 
time interval dt the following exchange events may take place. On the one hand, positive 
particles can exchange with adjacent holes to their right, negative particles with adjacent 
holes to their left, and positive particles with adjacent negative particles to their right, at 
rates 1, 1, and g, respectively. That is, between two adjacent sites one has 

+ 0^0+, 0-^-0, +-^-+, (1.1) 

with probabilities dt, dt and qdt. On the other hand, positive particles are introduced at 
the left end with rate a and leave the system at the right end with rate and conversely 
for negative particles. Thus, at site 1 

0^+, -^0, (1.2) 

with probabilities adt and pdt respectively; at site N 

0^-, +^0, (1.3) 
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with probabilities adt and pdt respectively. The dynamical rules of the model are symmet- 
ric under the simultaneous interchange of the sign of particles and of left and right, i.e., the 
dynamics of the positive particles moving to the right is the same as that of the negative 
particles moving to the left. Therefore, as long as this symmetry is not spontaneously 
broken, the densities and currents of the positive and negative particles are expected to 
be equal. 

The model has been studied in [4, 5] both in a mean field approximation and by sim- 
ulation. The mean field phase diagram exhibits a broken symmetry phase when the rate 
/? is smaller than a critical value Pc{c(,q)] in this phase the currents and the densities of 
the positive and negative particles are unequal. Numerical simulations of the model give 
similar results: the system evolves to one of two possible long lived states. These two states 
are related by the symmetry of simultaneous sign of particles and left-right interchange. 
In one of them the current and density of positive particles are larger than those of the 
negative particles, and in the other these inequalities are reversed. Outside of small regions 
near the boundaries, the states appear to be spatially uniform, and in the N ^ oo limit 
they become true steady states, hereafter also referred to as phases, to be identified with 
those of the infinite system studied, in the case = 1, in [6]. For any finite size N the sys- 
tem flips between these two possible states as it evolves, resulting in an overall symmetric 
state. This dynamics is characterized by a time scale (T), the average time between two 
consecutive flips, which diverges in the N oo limit. In [4, 5] it has been suggested, using 
numerical simulations and some simple considerations, that (T) grows exponentially with 
N. The flipping process itself takes place over a time interval much less than (T). 

To summarize, the two-species exclusion process studied in [4, 5] is known, at least on 
a heuristic basis, to exhibit the following two properties (when /3 < Pc)' 

(i) the system spends most of its time in one of the two symmetry related long lived states 
described above, each corresponding, in the N ^ oo limit, to a spatially homogeneous 
phase over the entire system, 

(ii) the average time between two consecutive flips, i.e., the average lifetime of a long lived 
state, is exponentially increasing with the size of the system: 

(T) ~ exp(iV/x). (1.4) 

The positive constant /j, in (1.4), hereafter referred to as the mass, is the analogue in the 
present situation of the barrier height, or activation energy, per site, characteristic of the 
law of Arrhenius for systems at equilibrium. 

These two properties, in analogy with the situation which prevails in the case of the 
two dimensional Ising model below the critical temperature, lead us to say that the model 



3 



exhibits spontaneous symmetry breaking. The analogy between the two models may be 
pursued. The escape rate j3 in the exclusion process plays the role of temperature, a and q 
being kept fixed. In particular, density fluctuations increase with and a transition from 
a broken symmetry phase to a symmetric phase is found when (3 = (3c- Finally changing 
the input and output rates at the boundaries to produce an asymmetry between positive 
and negative particles is akin to introducing a magnetic field in the Ising model (see the 
conclusion). 

The aim of the present paper is to demonstrate the existence of spontaneous symmetry 
breaking in the /3 ^ limit of the two-species exclusion model, by showing that the 
two properties mentioned above hold. In this limit the model can be mapped onto a 
much simpler toy model, for which the set of available configurations is restricted to those 
composed of three blocks, with j negative particles on the left, i/ holes in the middle and 
k positive particles on the right, with N = j + k + v (see eq. (2.1)). The dynamics of these 
configurations may be represented by a biased random walk in the positive quarter plane 
{j, k), with boundary conditions to be defined below. We demonstrate that the toy model 
fiips between two long lived states. These states coincide, except near the boundaries, 
with one or the other pure phase of the system, i.e., typically the larger of the (— ) and 
(+) blocks fills the system except for a region whose size does not grow with N: with high 
probability, max(j. A;) = N — 0(1). We are able to compute the average time between 
two consecutive flips, and flnd an exponential growth with the system size according to 
eq. (1.4). 

Let us finally emphasize that both properties: 

(i) existence of pure phases (related by symmetry), 

(ii) divergence of the time scale (T) with the system size, 

are necessary in order to claim spontaneous symmetry breaking, neither of them taken 
separately being sufficient. 

This may be illustrated by the following case. In the one-species exclusion process on an 
open chain, with input rate a and output rate studied in [7-10], the typical configuration 
on the phase coexistence line < a = P < 1/2 consists of two regions of size 0{N), one in 
each phase (one of low density a, the other of high density 1 — a), separated by a shock. 
The shock diffuses, hence its position is uniformly distributed over the chain in the steady 
state, and there is a characteristic time of order A^^ for the shock to traverse the system 
and hence for the system to change from one pure phase to another, without ever staying 
in one of them. In other words, there are no pure long lived states in this case; the steady 
state is a mixture. In conclusion, though the time scale (T) necessary for the system to 
pass from one phase to the other (without staying in one of them) diverges with the size. 
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no spontaneous symmetry breaking occurs in this case. Note that the N'^ divergence of 
(T) , instead of an exponential one, comes together with the absence of pure phases. This 
situation is presumably generic. We will come back to this example in section 5 since it 
corresponds to a limiting case of the two-species process. 

This paper is organized as follows. In section 2 we present simple physical arguments 
suggesting that, in the (3^0 limit, the dynamics of the exclusion process may be mapped 
onto that of the toy model. This correspondence is made rigorous in Appendix A. Section 3 
provides a description of the toy model as a biased random walk in the positive quarter 
plane. The equations for the stationary probability are given in section 4, as well as those 
for the average time between two flips. The latter involves two intermediate functions 
Uj^k and Vj^k, which have simple probabilistic interpretations. Appendix B provides a 
convenient way to compute all these quantities numerically. Section 5 is devoted to the 
study of two simple limiting cases of the random walk which are helpful to understand the 
general case. The analysis of the general case is given in section 6, where the asymptotic 
expressions of the stationary probability measure and of the average time between two 
flips are derived. We thus show that the steady state in the toy model is concentrated 
on conflgurations in which the system is almost entirely in a single phase. The analytic 
expression of the mass /i is given. In section 7 we determine the scaling behavior of Uj^k, 
Vj^k and the average time between two flips in the continuum limit, corresponding to a 
small bias and large distances. A summary and a discussion are given in section 8. 

2 Definition of the toy model 

In this section we show by means of simple physical arguments that in the limit where the 
rate /3 at which particles leave the system is small compared to the other rates defining the 
dynamics — a (particle input), q (particles interchange) and 1 (particle- hole exchange) — 
the exclusion model reviewed in the Introduction may be mapped onto a simpler model, 
which we term the toy model. A rigorous proof of the equivalence of the two models is 
provided in Appendix A. 

The time scale for the toy model is determined by the escape rate f3, and we therefore 
introduce a rescaled time r = 2(3t. In the (3^0 limit the only relevant configurations are 
those composed of three blocks, with negative particles to the left, holes in the middle and 
positive particles to the right: a typical configuration is of the form 
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All other configurations may be neglected. This can be understood as follows. 

First, since (3 is small, outputs of particles are rare on the time scale of the exclu- 
sion model. Any configuration different from those above will therefore rearrange itself — 
quickly, on the r time scale — until all the particles are waiting to exit. Suppose now that 
a particle exits. Either it will be replaced by a particle of the opposite sign, which will 
quickly traverse the system until it reaches the block of particles of its sign, or the hole 
thus created will itself travel through the system to join the block of holes in the middle. 
On the r time scale, then, the dynamics of the system involves transitions between various 
three-block configurations. These transitions may be described only in terms of the events 
taking place on the boundaries. Note that if the particle which exits is the last of its type, 
then the system will quickly fill with particles of the other type, leaving a configuration 
with only one block. Moreover, should the system empty completely, particles will quickly 
enter and rearrange to the form (2.1). Thus all configurations with j = or A; = 0, except 
for those with k = N or j = N, respectively, may be neglected in the /? ^ limit. 

Let us now give a quantitative description of the dynamics in this configuration space. 
We denote by (j, k) a three-block configuration, according to eq. (2.1). The integers j and 
k may take the values 0,1, . . . , N, with the restriction that j + k < N. Configurations 
(0, k) with <k < N and (j, 0) with < j < N do not occur. 

Consider now what may happen during a time interval At which is short on the r time 
scale but long on the t scale: P <^ At <^ 1. We can ignore the possibility that more 
than one particle exits in such an interval. If j > 1, then there is a small probability 
(1/2) At = PAt that the negative particle at site 1 leaves the system during this time 
interval. Suppose that this happens. 

(i) If J > 1, then immediately (on the t scale) either the hole thus created is filled by a 
positive particle, or the negative particle on site 2 exchanges position with the hole; the 
relative probabilities of these events are a/(l + a) and 1/(1 -|- a). Either event will be 
followed by a fast reordering of the particles, resulting in the configuration {j — 1, k + 1) 
or {j — 1, fc), respectively, unless j = N, when in the second case an additional negative 
particle will enter once the hole reaches the right end of the interval, resulting in the final 
configuration (A^, 0), that is, in no net change. 

(ii) If J = 1, then when the negative particle at site 1 leaves the system, the negative 
block disappears altogether. In this case positive particles are quickly introduced into the 
system at the left end and the system becomes filled with positive particles. The resulting 
configuration is (0, A^). 

Similar processes take place at the right end when a positive particle leaves the system. 
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To summarize these results we introduce rescaled rates 



which satisfy 

a + b = 1/2. (2.3) 
Then during the time interval dr the following stochastic changes take place in the system: 





(i-l,/c + l) 


ifj >1, 


with probability 


bdr, 


ij,k) - 


^ U-hk) 


if iV > j > 1, 


with probability 


adr, 




- (0,iV), 




with probability 


(1/2) dr, 


< 

i3,k) - 


^ U + i,k-i) 


if /e > 1, 


with probability 


bdr, 


ihk) - 


- (J,A:-1) 


if iV > A;> 1, 


with probability 


a dr, 




^ (iV,0), 




with probability 


(1/2) dr. 



(2.4) 



Let us emphasize again that states {j, 0) with j < N and (0, k) with k < N do not occur. 

These rules completely specify the toy model. They describe a random walk in the {j, k) 
plane, biased towards the south-west direction, that we study hereafter. 

3 The toy model as a random walk 

The rest of this paper is devoted to a quantitative analysis of the random walk defined 
by the dynamical rules (2.4) of the toy model. For the sake of clarity we summarize its 
definition first. 

A random walker hops from site to site on a two dimensional lattice, restricted to the 
positive quarter plane. The position of the walker is denoted by (j, k). During the time 
interval dr it hops from {j, k) to the south {j, k — 1) or to the west {j — l,k), each with 
probability adr. It may also hop to the north-west (j — 1, A; + 1) or to the south-east 
{j + l,k — 1) each with probability bdr {a + b = 1/2). When the walker reaches the j or 
the k axis it restarts instantly at site (AT, 0) or (0, AT), respectively. The random walker is 
thus subjected to non local boundary conditions. For short we will say that the axes are 
bouncing boundaries. 

Let us determine the velocity vector and the diflFusion tensor of this random walk. Con- 
sider the probability Pj,k{'r) to find the walker on site {j-,k) at time r. Its evolution in 
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time is given by the following master equation, obtained by conditioning on the last step 
of the walker 

^Pj,k{^) = a{Pj+i,k{^)+Pj,k+i{^)) + b{Pj+hk-i{^) + Pj-i,k+i{r)) -Pj,k{r), (3.1) 
leaving aside boundary conditions. This equation may be rewritten as 

^p(x, T) = ^a^ [p(x - ei, r) - p(x, r)] , (3.2) 

i 

with X = (j, /c), and where the (i = 1, . . .4) correspond to the four possible moves. To 
the south: ei — (0, —1), to the west: 62 = ( — 1, 0), to the north-west: 63 = (—1, 1), to the 
south-east: 64 = (1, —1), with rates ai = 02 = a, and as = 04 = 6 respectively. 

The velocity vector reads [11, 12] 

V= lim -(x(t)) = y a^e^, (3.3) 

r— »oo T 

i 

while the components of the diffusion tensor D are given by 

i 

In the present case one finds 

V =(:::), (3.5) 

and 

^ 1 / 1 - a 2a - 1 \ 

The eigenvalues of this matrix yield the two diffusion coefficients 

D, = ^, D„=^ (3.7) 

where perpendicular and parallel refer to the south-west direction of the velocity (or bias). 
Note that the velocity and the parallel diffusion coefficient vanish as a ^ 0. 

The toy model captures the essence of the original exclusion process, when (3 is small; 
in particular the mechanism leading to spontaneous symmetry breaking has an intuitive 
explanation in terms of the random walk. Let us first review this mechanism for the 
exclusion process. Suppose the system is in the state of high density of (+) particles. In 
this phase there is a low flux of negative particles and holes entering the system at the 



8 



right end and leaving at the left end. As the system evolves, a block of negative particles 
may temporarily be formed at the left end due to some fluctuation. As a consequence holes 
entering the system become trapped between the positive and negative regions. Usually, 
the block at the left end leaves the system after some time and the system relaxes to the 
all (+) state. However, if the block persists long enough, the system may be filled with 
holes, and thus has a chance of flipping to the all (— ) state. 

This mechanism has also a clear interpretation in the framework of the biased random 
walk of the toy model. Consider the walker starting at (0, A^), i.e. in the all (+) state. 
With large probability the walker will not go far east against the bias. It will therefore soon 
hit the k axis and bounce to its starting point (0, N). However, with a small probability, 
vanishing as A ^ oo, it may reach the other end of the system (A, 0). It is thus clear that 
the walker spends most of its time moving near one of the two endpoints of the system, 
occasionally traversing to the other end. Therefore the stationary probability measure of 
the walker is concentrated around both endpoints (0, A) and (A, 0). The average time for 
the walker to move from one end of the system to the other is the average time between 
two flips in the toy model. This time is denoted hereafter by (Tq). It is related to the 
average time between two flips of the original exclusion process in the /? — > limit by 

(T)^^^ (3.8) 

The stationary probability measure as well as the average time between two flips in the 
toy model are studied in the next sections. 



4 Equations for the stationary probability and for the average time between 
two flips 

Two questions are to be answered: 

(i) What is the probability pj^k for the walker to be at site {j, k) in the steady state (i.e. 
the stationary probability measure)? 

(ii) What is the average time Dj^^ for the walker, starting from site (j, /c), to reach any 
site on the j axis for the first time? The average time between two fiips in the toy model 
is simply 

(To) = L»o,iv. (4.1) 

In this section we establish the equations fulfilled by these two quantities. Their so- 
lutions, mostly in the asymptotic A — > oo regime of interest, are presented in the next 
sections. 



9 



4-1 The stationary probability measure 

In the steady state, the equation for the probabihty pj^k to find the walker on site {j, k) 

is 

Pj,k = a iPj,k+i +Pj+i,k) + b (pj+i^k-i +Pj-i,k+i) (4.2) 
as a simple consequence of eq. (3.1), with boundary conditions 

p,- = (j^ N) 

Po,k = (k^N) (4.3) 
P3,k = ij + k>N). 

The first two conditions are the only relevant consequences of the bouncing property of 
the axes. 

The stationary probability measure obeys the symmetry property 

Pj,k^Pk,j- (4.4) 
In addition the following normalization condition is imposed 

E^'^-'^ = 1- (4.5) 

j,k 

4-2 The average time between two flips 

Dj k, as defined above, is the first-passage time through the j axis. It is therefore 
convenient to consider this axis as absorbing (the k axis is still bouncing). Conditioning 
on the first step of the walker, during the infinitesimal time interval dr, yields 

Dj^k = adT{Dj^k-i + dr) + adT{Dj-i^k + dr) 

+ bdT{Dj+i^k-i + dr) + bdT(Dj_i^k+i + dr) (4.6) 
+ [l-{2a + 2b)dT]{Dj,k + dT), 

hence 

Dj^k = 1 + a{Dj^k-i + Dj-i,k) + b{Dj+i,k-i + Dj-i,fc+i), (4.7) 
where j and k = 1, . . .,N. The boundary conditions read 

Do,k = Do,N (4.8a) 
Dj,o = (4.8b) 
bDo,N = bD^^N-i + 1. (4.8c) 
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The first equation reflects the non local bouncing property of the k axis, while the second 
one expresses that the j axis is absorbing. The last equation comes from the fact that, 
starting from site (0, A?"), the walker has probability bdr to hop to site {1,N —1), and the 
complementary probability 1 — 6 dr to stay on this site. Note that the recursion relations 
found for the -Dj,fcj eqs. (4.7), are the same as those for a discrete-time process [13], since 
the hopping rates of the walker add to 1. 

In order to unravel the cumbersome condition eq. (4.8a) imposed by the bouncing bound- 
ary, it is convenient to introduce two functions uj^k and vj^k satisfying 

Uj^k = 1 + a{uj^k-i + Uj-i^k) + b{uj+i,k-i + Wj-i,fe+i) (4.9a) 
Vj,k = divj^k-i + iJj-i,k) + b{v3+i,k-i + Vj-i^k+i), (4.9b) 

with the boundary conditions 

UQ,k = 0, vo,fc = 1 (4.10a) 
Uj,o = 0, Vj,o = 0. (4.10b) 

Uj^k and Vj^k obey the symmetry properties 

Uj,k = Uk,j, Vj^k + Vk,j = 1. (4.11) 
Then the solution to eqs. (4.7) and (4.8a-c) is given by 

Dj,k = Uj^k + Dq^n Vj^k, (4.12) 

and in particular by 

D„,« = I^Mtl+V^. ,4.13) 

1 — 1^1,iV-l 

We are actually interested in determining Do, at in the asymptotic regime, physically 
relevant, where N is large. As shown by eq. (4.13) two independent sets of equations for 
Uj^k and Vj^k have to be solved separately. This problem is studied in the next sections. 

We conclude by giving simple interpretations to the quantities met in this section and 
to the relations between them, in particular to the key eq. (4.13). 

(i) As can be deduced by considering the equations fulfilled by Uj^ki this quantity represents 
the average time for the walker starting from (j, k) to reach either of the axes for the first 
time. Similarly Vj^k is the probability for the walker starting from {j, k) to reach the k 
axis without hitting the j axis. Generalizing the language of the classical Gambler's ruin 
problem, Uj^k is the duration of the game, and Vj^k the probability of ruin: the walker is 
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ruined when reaching the k axis, whereas it wins when reaching the j axis. The game stops 
in both cases. When considering the quantities Uj^^ and Vj^j- it is therefore convenient to 
consider both axes as absorbing. 

Hence equation (4.13) also gets a simple interpretation. Its denominator is the prob- 
ability for the walker, starting from (1, N — 1) to reach the j axis, that we may name 
the probability of success. In the numerator, 1/b is the average time for the walker start- 
ing from (0, A^) to reach site (1, A'" — 1), while wi,Ar-i is the average time, starting from 
(1, — 1), to reach either of the axes. Altogether the numerator of (4.13) represents the 
average duration of an attempt for the walker starting from (0, A^): once the walker has 
reached site (1, A^ — 1), either it succeeds with probability 1 — vi^n-i in traversing the 
quarter plane, or it fails with probability vi^n-i, returning to its starting point (see also 
eq. (4.19) below). Alternatively this duration represents the average time between two 
bounces on either of the axes. Denoting this duration by (7b), we have 



and (4.13) reads, using (4.1) 



in) = ^ + ui,N-i, (4.14) 



(To) = ^ ^^^^ . (4.15) 

1 - Viw-l 



The probability of success is therefore the ratio of the frequency of flips (or success) to the 
frequency of bounces (or attempts) . 

(ii) The average time between two bounces (Tb) can be independently evaluated as follows. 
Just after a bounce the walker either passes through the bond from (0, A^) to (1, A^ — 1), 
or through the bond from (A^, 0) to (A^ — 1, 1). Hence (Tb) = 1/(2J), where J — 6po,Ar is 
the probability current through either of the two bonds. Thus 

(^^b) = (4.16) 



Comparing eqs. (4.14) and (4.16) yields 

1 



which will be used hereafter. 



= l + 6«i,iv-i (4.17) 



(iii) Finally, an alternative interpretation of (4.15) is provided by the following considera- 
tions. The average number of attempts of the walker until the first success is 

(m) = V m ivi,N-iril - vi,N-i) = , ""''^"^ • (4.18) 

^0 l-^l,iV-l 
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Denote by u^i\_i the average time for a walker starting from {1,N — 1) to touch the j 
axis without 
and k. Then 



axis without touching the k defined similarly, by exchanging the roles of j 



(To) = (m) Q + u?],_,^ + Q + , (4.19) 

where the first parenthesis represents the average duration of a failure, the walker starting 
from (0, A^) and touching the k axis first, while the second one represents the average 
duration of a fiip (or success), the walker starting from (0, N) and touching the j axis first. 
Combining eqs. (4.18, 4.19) and noting that 



(Tb) = - + u,,N-i = ^i,Ar-i (^^ + «i:Ar-i j + (1 - Vi,N-i) [ ^ + <jv-i ) (4.20) 
again yields the expression (4.15) of (To). 



5 Two limiting cases 



In this section we consider two simple limiting cases a = and 6 = 0, which allow 
for exact closed-form solutions, that will be helpful to understand the general case. We 
use here the discrete-time language, for simplicity, since it yields the same equations and 
results as the continuous-time formalism used so far, because the hopping rates of the 
walker add to 1. 



5.1 The case a = 

In this limiting case the toy model reduces to the classical Gambler's ruin problem on 
the segment j + k = N , with equal probabilities 6=1/2 for the walker to hop to the right 
or to the left, and with refiections on sites {0, N) and (A'", 0). The stationary probability 
measure is uniform on that segment, i.e.. 

On the other hand, eqs. (4.10a) and (4.10b) have the following simple solutions 

k 

Uj,k^jk, Vj,k = j—j^- (5.2) 
By inserting these results into eqs. (4.13) and (4.14), we are left with 

(Tb) = N + 1, (To) = Do,N = NiN + 1). (5.3) 



13 



The asymptotic N'^ law is due to the diffusive nature of the process. One also finds 



Dj,k = N{N+l)-j{j + l). (5.4) 

5.2 The case b = 

This case also yields a classical problem, Le probleme des points, considered by Fermat 
and Pascal [14]. 

In this case, once the random walker has reached either of the axes, it stays there. There- 
fore for 6 = the model is not ergodic even for a finite system size, and any probability 
measure supported on the two points (0, A^) and (A^, 0), i.e., satisfying po,iv +Pjv,o = 1) 
is stationary. In the 6^0 limit of the toy model we obtain the symmetric measure, i.e., 
Po,N = Pn,o = 1/2. 

Let us now compute Uj^k and Vj^k- Consider the more general case where the probability 
of hopping south is p and of hopping west is r, with p + r = 1. Eqs. (4.10) and (4.11) can 
be solved by means of generating functions. Dealing first with the vj^k, we introduce the 
two-variable generating series 

CO oo 

V(x,|/) = J] x^' (5.5) 

j=i fe=i 

Eq. (4.10b) yields a closed equation for the function V{x, y), which can be solved explicitly. 
We thus obtain the rational form 

V{x, y) = (5.6) 

{)- - y){)- - rx - py) 

By expanding this expression as a double series, we obtain 

k-i 



n 

n=0 



Note that the nth term in the sum represents the probability that the walker, starting at 
(j, /c), reaches the k axis after having made n < k steps south. In particular one gets 

vi,fc = l-/, Vj,i=r'. (5.8) 

In a similar way, eq. (4.10a) yields the following rational expression for the generating 
series of the Uj^k 

U(x,y) = -, r- r. (5.9) 

^ '^^ {l-x)il-y){l-rx-py) ^ ^ 
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By expanding this result we get 
and especially 

ui,k = , = • (5.11) 

r p 

(Tb) and (Tq) are divergent when 6^0. Indeed by inserting the above results for 
p = r = 1/2 into eqs. (4.13), (4.14), we obtain their asymptotic behavior in the 6 — > 
regime: 

(Tb)^^, (To)^^^ (6-0). (5.12) 

The second result has the form of eq. (1.4), with a mass = ln2. 

Coming back to the original two-species exclusion process, one notes that the two cases 
considered here correspond to two extreme situations. 

The case a = corresponds to a infinite. In this case the stationary state is known [5] . 
No holes are present in the system and the remaining (+) and (— ) particles play the role of 
particles and holes, respectively, in the equivalent one-species totally asymmetric exclusion 
process to which the original process reduces. Moreover, in this equivalent process, the 
rates for input and output of particles are equal to a = P/q, implying — /3 being small — 
that the system is at coexistence between a low density phase with a density of particles 
equal to a and a high density phase with a density of particles equal to 1 — a, as reviewed 
in the Introduction. The instantaneous configuration of the system consists of two regions 
of order A^, one of low density of (+) (or high density of (— )) on the left, separated from 
another one of high density of (+) (or low density of (— )) on the right by a shock. The 
location of the shock diffuses between and A^. In the toy model, the location of the 
random walker is precisely that of the shock. It is uniformly distributed over the whole 
system. As said in the Introduction this case does not correspond to spontaneous symmetry 
breaking. 

In the other extreme case 6=0, i.e., a = 0, the dynamics stops in one of the states 
of broken symmetry. Hence this situation does not correspond to genuine spontaneous 
symmetry breaking. 

6 The general case 

We now turn to the analysis of the large- A?" behavior of the stationary probabilities pj^k 



15 



and of the average time (Tq) between two flips in the general situation (0 < a < 1/2, 
6= 1/2 -a). 

6.1 The stationary probability measure 

The stationary probabihty measure pj^k is the solution of eq. (4.2), with the boundary 
conditions (4.3). It is convenient to use as co-ordinates j and the number of holes u, such 
that j + k + u = N. With the notation p^^'^^ = pj^^, eq. (4.2) reads 

pf^ = a (p^-'^ + p%-'^) + b (p^i + p^%) . (6.1) 

with boundary conditions inherited from eq. (4.3). It is shown in Appendix B how this 
equation can be solved numerically in a recursive way, according to increasing values of u, 
starting from u = 0. 

In the asymptotic N ^ oo regime of interest, the probability measure is expected to be 
concentrated in two finite regions around the points (0, A^) and (A^, 0), each region being 
decoupled from the other one. This phenomenon, already observed in section 5 for 6 = 0, 
is a consequence of the bias present in the dynamical rules of the toy model. 

Because of the symmetry property (4.4), it is sufficient to consider the vicinity of the 
point (0, A^). We are thus led to look for an asymptotic solution p^^^ of eq. (6.1), inde- 
pendent of N in the A" — > oo limit, and localized, namely decreasing as either j or i/ gets 
large, i.e., for 1 <^ j, <^ A". The normalization of the probability measure now reads 

i/>0 j>0 

since the two decoupled regions around (0, A^) and around (A^, 0) bear equal weights. 

In order to study these asymptotic probabilities Pj'^\ it is convenient to introduce the 
one-variable generating functions 

Eq. (6.1) then leads to the recursion relation 

Pj{x) = ax[Pj{x) + Pj+i{x)] +b[Pj+i{x) + Pj-i{x)], (6.4) 
with (see eq. (4.3)) 

Po{x)=pi^l (6.5) 
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The general solution of eq. (6.4) reads 

Pj{x) = A{x)[z_{x)Y + B{x)[z+{x)Y, (6.6) 
where z± (x) are the two solutions of the characteristic equation 

{ax + b)z'^ - {l-ax)z + b = 0, (6.7) 

hence 

, , 1 — ax ± \/ A(x) ,„ 

^±(^) = — 7m ' (^-^^ 

l[ax + b) 

where 

A(a;) = a^x^ - 4a(l - a)x + 4a(l - a). (6.9) 
This expression vanishes for a; = a;i and x = x^-, with 

xi = ^ (^1 - a - V(l-a)(l-2a)) , X2 = ^ (^1 - a + ^(1 - a)(l - 2a)) . (6.10) 
The associated values of 2; read 

1 - 2a + 2^(1 -a)(l -2a) 



2:1 = 2;±(a;i) = 

Z2 = Z±{X2) = 



3 -2a 

1 - 2a - 2y/{l- a){l-2a) 
3 -2a ' 



(6.11) 



The branch Z-{x) is monotonically increasing from Z- = exp(— u) for a; = 0, with the 
notation 

1 1 . /, X 1 ± x/4a(l - a) 

cosher = — - = , I.e., expiicr) = (6.12) 

26 l-2a' ^ ' l-2a ^ ' 

to Z- = Zl for X = xi. Conversely, z+{x) is monotonically decreasing from z+ — exp(a") 
for X = to z+ = for x = xi. 

Since the generating functions Pj{x) are analytic at least in a neighborhood of the origin 
X — 0, and decrease to zero for large j, the second term in the general solution (6.6) is 
ruled out. Taking into account the boundary condition (6.5) at j = 0, we are left with 

P,{x)=p^^\z_{x)y. (6.13) 

Then, setting a; = 1 in eq. (6.13), and using Z-{1) = 1 — 2a, yields 

i/>0 
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By summing this expression over j, and using eq. (6.2), one gets the very simple result 

p'^Q^ = a. (6.15) 

The average time between two bounces (Tb) in the N ^ oo regime can then be evaluated 
by means of eq. (4.16). We thus obtain 

(T.) = ^. (6.16) 

Eq. (6.13) is the main result of this section, from which the following consequences can 
be derived. 

(i) A quantitative characterization of the concentration of the probability measure around 
the point (0, N) (and similarly around {N, 0)) is provided by the following order parameter 

M2(iV) = <(^)'>, (6.17) 



defined as an average with respect to the stationary measure. The existence of the asymp- 



totic probabilities p[^\ independent of N, implies 



M2{N) = l-^ + o(^], withA = 4(j) + 2(zy), (6.18) 



the averages being taken with respect to the Pj'^\ The result (6.13) yields 

(j) = 2^jP,(l) = ^^ (6.19) 



and 



j>0 V / x=l « 

so that we have 

^^2(l-2a)(2-a) 

a 

Note the simple interpretation of (6.19): (j), the average number of (— ), i.e., of minority 
particles in the steady state, is equal to the input rate a. It is also equal to Dj;y/\4, using 
eqs. (3.5, 3.6). 

(ii) From the expression (6.13) of Pj{x) it is also possible to extract the p^^^^ by means of 
the following contour integral encircling the origin in the complex x-plane 
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We will restrict ourselves to the asymptotic behavior of these probabilities, when both j 
and v are large {1 <^ j^v <^ N), at a fixed angle < (p < tt/A between the negative A;-axis 
and the vector joining (0, N) to {j, k), namely we set 

N — k — j + i' = p cos 09, 

(6.23) 

j = p sin if. 
We are thus led to evaluate the integral 

— ^ ^ 2^~~ ^"^P { ~ ^^^^^ ~ f^) In a; — sin (p In z- (a;)] | , (6.24) 

for (f fixed, in the p — > oo limit. The direction ip = corresponds to — > oo at constant 
j, whereas ip = tt/A corresponds to j — >^ oo at constant v. 

A saddle-point calculation leads to the estimate 

p"^^ ~ exp {-pp{ip)} , (6.25) 

with 

= (cos<^ — sin (^) In — sin In 2;_(a;c)) (6.26) 
and the saddle point Xc{f) is given by 

We have thus demonstrated that the asymptotic probabilities pj^'' decay exponentially in 
all directions away from the point (0, A^), with an inverse decay length given by the mass 
p,{(f) in the direction defined by the angle (p. Eqs. (6.26) and (6.27) provide a parametric 
representation of p,{ip). It increases monotonically between <^ = 0, where Xc = xi, and 

^(0) = = Inxi, (6.28) 

with the notation (6.10), and p — 7r/4, where Xc — 0, and p{n/4) — a"/\/2, with the 
notation (6.12). 



6.2 The average time between two flips (Tq) 

This quantity is given by eq. (4.13). First, the asymptotic behavior of wi^jv-i for large 
N is easy to determine. Indeed the walker will visit the j-axis with a vanishingly small 
probability. One may therefore neglect the presence of this boundary, and just consider 
the displacements of the walker parallel to the j-axis. Thus Uj^k reduces to uj, the average 
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time for the walker starting at site j to reach the boundary, i.e., site j = 0. During dr the 
walker may either hop to the left with probability (a + 6)dT or to the right with probability 
bdr or stay on site j with probability adr. Therefore the set of Uj obeys the equation 

Uj = 1 + {a + h)uj-i + huj j^i + auj, (6.29) 

with the boundary condition uq = 0. Equivalently eq. (6.29) is obtained by dropping the 
^-dependency in eq. (4.10a). The solution to eq. (6.29) is simply Uj = j/a. Coming back 
to the original two-dimensional situation, we find 

«i,iv-i^- (iV>l). (6.30) 
d 

This result provides an alternative derivation of eq. (6.15), by means of eq. (4.17). 

The analysis of the asymptotic behavior of i is more difficult. Introducing the gener- 
ating function V{x,y) as in eq. (5.5), we see that eq. (4.10b) is equivalent to 

2 2 

D{x,y)V{x,y) + bxy[xVi{x)+yV2{y)] = ^^^^ (6-31) 



where 



X-' , 



V2{y) = ^{x,y)a;=o = Xl^^i.fe/, 

fc=l 



(6.32) 



and 



D{x, y) — xy — a{x + y)xy — b{x'^ + y^). (6.33) 



The characteristic equation D{x,y) = defines the characteristic curve of the problem. 
It is a rational cubic curve, with the origin as a double point, and the three lines x = —b/a, 
y = —b/a, and x + y = {1 + 2b)/a as asymptotes. Figure 1 depicts the characteristic curve 
for a = 0.4, 6 = 0.5- a = 0.1. 

We notice that the characteristic equation is equivalent to eq. (6.7), up to the substitution 
y = xz. Solving the characteristic equation for y, we thus get 

y±{x) ^ xz±{x), (6.34) 

with the notation (6.8). 

The characteristic curve has either a horizontal or a vertical tangent at the four points 
marked on Figure 1. Their co-ordinates read A = {xi,yi), B = {yi,xi), C = {x2,y2), 
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D = {y2,X2), where j/i = xiZi, j/2 = X2Z2, with the notations (6.10) and (6.11). Figure 
1 also provides a geometric illustration of the saddle-point approach to the stationary- 
probability, presented in section 6.1. Indeed, as the angle (f varies from to it/ 4, the 
saddle point runs along the dispersion curve, in co-ordinates (xc, j/c) = {xc,XcZ-{xc)), 
from the point A to the origin O. 

Let us come back to the Vj^k- Since these numbers are bounded, the function V{x, y) is 
analytic in a complex domain V containing at least the product of disks < 1, \y\ < 1. 
Whenever x and y are in V, and related by the characteristic equation D{x, y) = 0, 
eq. (6.31) simplifies to 

^Vii^) + y^M = ^^^^y (6.35) 

This equation shows that at least one of the two functions Vi{x) or V2{y) becomes singular 
when the relation between x and y becomes itself singular, i.e., at the marked points 
mentioned above, as may be seen by taking the derivatives of both sides. 

The points A and B are the relevant ones, since they are the closest to the origin. Since 
these two points play symmetric roles, we can consider point A for definiteness. We have 
j/i < 1 < xi, so that Vi{x) is singular at a; = a;i, whereas V2{y) is regular at y = j/i. 

We set 

x = xi-€x, y = yi + €y, (6.36) 
so that the characteristic equation reads, to leading order in and 

e^^(^^l±^) i, with w = J(l-a)(l-2a). (6.37) 
\ a;i y aw 

By inserting this estimate into eq. (6.35), and expanding both sides to leading order in e^, 
we are led to the conclusion that Vi (x) has a singular part of the form 



Vi{x) « Vi{xi) - K^xi - X, (6.38) 

with 



d 



ax\ + bdy 



yV2{y) - 



26(1 -y) 



(6.39) 



y=yi 



By inserting the singular behavior (6.38) into the definition (6.32) of Vi{x), we get the 
estimate 

vj,i ^ ^^^pW2) j-3/2 exp(-^7,) » 1), (6.40) 
where /i has been defined in eq. (6.28). 
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We can now insert both estimates (6.30) and (6.40) into the expression (4.13) of (Tq). 
We thus obtain 

(To) ^ L N^/^ ex.p{Nii) (N > 1), (6.41) 



with 



L = 



20F 



(6.42) 



a(l - 2a)Kexp{3ii/2)' 



We have thus derived the full asymptotic expression of the average time between two flips. 

• The mass ^ is given by eqs. (6.10), (6.28). It vanishes for small values of a as |U = 
a/4 + llaV32 + ■ ■ • This regime corresponds to a continuum limit, to be investigated 
in detail in section 7. The opposite limit a ^ 1/2 yields fi In 2, in agreement with 
eq. (5.12). The mass increases monotonically between these two limits. For instance in 
the case a = 6 = 1/4 one finds = In (2(3 - a/6)) = 0.096 237. 

• The exponent 3/2 is universal, i.e. independent of a, except for the two limiting cases 
a = and b = 0, investigated in section 5. The general analysis presented above fails 
in these cases, because the characteristic curve is degenerate. We have indeed D{x,y) = 
— {x — y)^/2 for a = 0, and D{x,y) = xy{2 — x — y)/2 ior b = 0. 

• The absolute prefactor L only depends on a. It cannot be determined within the present 
approach. The analysis of the continuum limit will provide a quantitative estimate of the 
behavior of this prefactor as a — > 0. Indeed eq. (7.45) yields L f» 2/V'7T^a, i.e., K 'n'^ j 

7 The continuum limit 

In this section we determine the scaling behavior of Uj^k, vj^k and (Tq) in the continuum 
limit, corresponding to the regime of a small bias (a — > 0). It will turn out that this analysis 
provides a complete description of the crossover between the diffusive power law (5.3) and 
the exponential law (6.41). 

We consider first the probability of ruin Vj^k- We assume that Vj^k is a smooth function 
of the continuous rescaled variables X = ja and Y = ka. Expanding eq. (4.10b) to the 
lowest non-trivial order in a leads to 



dX 



d 



_d_ 
dY 



.(x.y) = 2(A + ^)„(x,n 



(7.1) 



hence to 



d^^ dd 



(7.2) 
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with 

'^^X-Y = {j- k)a, 

(7.3) 

e = X + Y = {j + k)a. 

We thus obtain a heat equation, or diffusion equation, where ^ and 9 are respectively the 
space- hke and the time-hke variables. Note that the time-like axis is oriented along the 
direction of the bias. This partial differential equation inherits its boundary conditions 
from eqs. (4.11a, b) 

^(^ = -^,^) = i, v{i = e,e) = Q. (7.4) 

These moving boundary conditions are diflBcult to handle. It is convenient to introduce 
the dimensionless space-like variable 

such that 

£, = {l-2x)e. (7.6) 

Eq. (7.2) thus transforms into 

0_,a-2.).|-4.^|.o. (7.7) 

with the boundary conditions 

^;(a; = 0, ^) = 0, i;(a; = 1,^) = 1. (7.8) 

The Gaussian 

G{x, 9) = (47r^)- V2 exp {-9{x- l/2f) (7.9) 

is an elementary solution of eq. (7.7), neglecting boundary conditions. We are therefore 
led to the change of function 

v{x, 9) = G{x, 9)Fy{x, 9), (7.10) 
which inserted into eq. (7.7) gives 

A last change on the time-like variable 
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leads to a heat equation for Fy , 
with boundary conditions 

\ 1/2 /I 



F,(a; = 0,£) = 0, = 1, e) = /^(e) = J exp (^-— J . (7.14) 

and with hme_»_oo Fy{x,e) = 0. 

We thus have to solve the heat equation on the unit interval, with prescribed time- 
dependent boundary conditions. The solution is given by the convolution 

nOO 

F,{x,e)= / ds' Ue-e')g,{x,s'), (7.15) 
Jo 

where 

gy{x, e) = 2n J^(-l)^"^ k sm{kTTx) exp{-k'^TT'^e) (7.16) 
fc=i 

is the appropriate Green's function of the heat equation (7.13). We mention for further 
reference its alternative expression, obtained from eq. (7.16) by the Poisson formula 

g,{x,e) = {47Te^)-^/^ (2m + 1 - x) exp (- ^^"^ + ^ ~ V (7.17) 

vn= — oo ^ ' 



We are thus left with the following explicit formula for the scaling behavior of the Vj^u 
in the continuum limit 

/"°° d£ f 6 \ 

v,,,^v{x,e) = exv{-e{x-l/2f) (i + 4g,)i/2 ^"P( 4(i + 4g,) j 9vix,'). (7.18) 

We recall that the scaling variables x and 9 are related to j and k by eqs. (7.3, 7.5). 

We are especially interested in the quantity vn-i,i, which enters the expression (4.13) 
of (To). Eq. (7.18) yields 

a^9) (d = Na), (7.19) 
where ^{0) is the following scaling function 

4W = -_.(. = 0,.) = -y^ ^^-^e.p(-_),„(e), (7.20) 
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with 

9o{£) = Q^9v{x = 0,£) 



= expl-pTT^e) 21) 

fc=i 

3^-l/2 /(2m + 1)2 \ / (2m + l)2\ 



The scahng behavior of the duration of the game Uj^k can be investigated in a similar 
way. Expanding eq. (4.10a) now leads to an inhomogeneous heat equation, 

with the moving boundary conditions 

u{(, = -e,e) = u{i = e,e) = Q. (7.23) 

We again use the dimensionless space-like variable setting 

u{x, e) = G{x, e)F^{x, e), (7.24) 

and we change the time-like variable from 9 to e. We thus obtain an inhomogeneous heat 
equation for F„, 

,7.25) 

with boundary conditions 

F^{x = 0, £) = F^{x = 1, £) = 0, (7.26) 
and with lim^^-oo Fu{x, s) = 0. The source term reads 

We now have to solve the heat equation on an interval, with a prescribed time-dependent 
source, and Dirichlet boundary conditions. The solution is again given by a convolution, 

Fu{x,e)= de' dyfu{y,s-e')gu{x,y,e'), (7.28) 
JO ^0 
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where 

oo 

gu{x, y,s) = 2 sin(A;7ra;) sin(A;7rj/) exp(— A;^7r^£) (7.29) 
fc=i 

is the appropriate Green's function of the heat equation. 

We are thus left with the following explicit formula for the scaling behavior of the uj^k 
in the continuum limit 



Uj 



2^2 

,fc ^ u{x, e) = —exp{- e{x - 1/2)2) 

(7.30) 



X 



L^'L (l + 4^e)V2-H4(rT4^J 



We are especially interested in the quantity W7v-i,i, which enters the expression (4.13) 
of (To). Eq. (7.30) yields 

«Ar-i,i -*(^) {9^Na), (7.31) 

ft 

where '^{9) is the following scaling function 

. ^ JL^. . 0..) = ^^.p (M^^ 

(7.32) 

where we have used the identity 

d 

—gu{x = 0,y,s) = gy{l-y,s). (7.33) 

Putting together both results (7.19, 7.31) we obtain the following expression for the 
scaling behavior of the average time between two flips 

(To)«^A(^) (e^Na), (7.34) 



a2 



where the scaling function A(^) reads 

MO) = (7.35) 

The regimes of small and large values of the time-like scaling variable 6 = Na are of 
special interest. We now show how the scaling law (7.34) describes the crossover between 
the diffusive law (5.3) in the unbiased case (^ <^ 1), and the exponential law (6.41) in the 
generic biased case {9 ^ 1). 

• ^< 1. 
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In this regime it is most convenient to consider eqs. (7.2, 7.22) and to look for solutions 
to them in the form of power series in d. We thus obtain 



u{x, e) = -^x{l - x) [1 - d + —{4: + X - x^) + ■ 

v{x, 6*) = a; + a;(l - a;)(l - 2a;) ( -- + — (1 - 2a; + 2a;^) + 

V 3 15 



(7.36) 



The leading terms of the above expressions agree with the results (5.2) corresponding to 
a = 0. The above results also yield the following expansions for the scaling functions 

i-e + — + ---y (7.37) 

The last of these expansions yields more explicitly 

(T„)«iV^(l-Hf5 + !I^ + ...) (N»l,Na«l). (7.38) 

The leading term reproduces the asymptotic behavior of the exact result (5.3), character- 
istic of diffusive motion. Note that the first correction term to this law is negative. This 
may be understood as follows. In the current regime (i.e., a » A?" ^ 1), we typically 
expect a small number of zeros to separate the (+) and (— ) blocks, thus leading to a small 
negative correction to the average time between two flips of the strictly diffusive shock 
behavior. 

• ^> 1. 

This regime is the most relevant physically, since it corresponds to the asymptotic be- 
havior (6.41). The asymptotic behavior of the scaling functions can be determined from 
the integral expressions (7.20, 7.32). 

The calculation for ^{9) is as follows. The double integral in eq. (7.32) is dominated by 
the region where both variables e and y are small. In this regime the function gy{l—y, e) can 
be approximated by the m — term in the sum (7.17). Choosing the rescaled integration 
variables z = y/ {9e) and C = 9^e, we are left with 
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The integral can be performed by changing variable from ( to rj = + ^C/O). After 
some algebra we are left with the very simple limit behavior 



*(^) ^1 > 1), (7.40) 

in agreement with the more general result (6.30). The correction to this limit vanishes 
exponentially as exp(— 0/4), since it is due to the upper bound ?7max = 6*74 of the 77-intcgral. 

Let us now compute $(0). We perform the change of variable ^ = 0/(1 + AOs) in the 
integral representation (7.20), with the last expression of eq. (7.21) for the function qq, 
obtaining thus 

2^cp(-0/4) 

*W ~ (7^03)1/2 ^' (7-41) 

with 

^=E / dC(2(2m + l)2C-l)exp[-((2m+l)2-l/4)C] 

^ (2m + l)^ + l/4 _.y-|^ 1 , 1 \ tt^ 
^o((2m+l)2-l/4)' ^oV(4™ + l)' {^rn + 2>Y) 4" 

We thus obtain the estimate 

^(^^ ^ \ {0 » 1). (7.43) 

The asymptotic estimates (7.40) and (7.43) yield 

3 /2 

A(0) ^ 2 exp(0/4) (0 > 1), (7.44) 



namely 



iV3\l/2 

(To) 2 ( — exp(iVa/4) {Na > 1). (7.45) 



aTT^ 



This result exhibits an exponential growth with mass = a/4 and a power-law prefactor 
in A?"^/^. These results are in agreement with the analysis of section 6. 

In conclusion we see that the analysis contained in this section fully describes the 
crossover from the diffusive behavior of a shock (0 <^ 1, i.e., a » A » 1) to the symmetry 
breaking regime (0^1, i.e., A ^ o; ^ 1). This crossover behavior is illustrated on Figure 
2, showing a logarithmic plot of (Tq) against 6 = Na. The numerical data for a = 0.1, 0.06, 
and 0.03, obtained by means of the approach described in Appendix B, smoothly converge 
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toward the result (7.35) of the continuum hmit. The smoothness of the corrections to the 
continuum hmit has aheady been noticed below eq. (6.42), for the mass /i. 

It is striking to observe the richness of the behavior of the original exclusion process: 
even in the regime studied here, i.e., /? — > and a — > oo, the curve obtained for (Tq) is 
very similar to that obtained in [4] by numerical simulations of the exclusion process with 
finite values of a and /? (see Fig. 4 of that work) . 

Finally we notice that the stationary probability measure pj^k also exhibits a scaling 
behavior in the continuum limit. A detailed investigation of this behavior is not needed. 
We just mention that eq. (4.17) implies the following scaling law for the probability at 
both endpoints (0, N) and {N, 0), in the regime a <^ 1, N 1 

Po,n^pT^^^ (e^Na). (7.46) 

As expected, this result interpolates between the a = limit, eq. (5.1), and the N ^ oo 
regime, eq. (6.15). 

8 Conclusion 

The analysis presented in this paper confirms the intuitive picture, given at the end of 
section 3, of the mechanism leading to spontaneous symmetry breaking in the toy model, 
and quantifies it in several respects. 

First, the stationary probability measure pj^k is localized in a symmetric way around 
the endpoints (0, N) and (iV, 0), in the asymptotic regime of a large system size N. This 
is intuitively clear: the bias of the random walker is in the south-west direction so that its 
drift is towards whichever boundary is nearest. It then becomes localized near one of the 
endpoints due to the bouncing property of the boundary. This localization is quantified 
by the behavior (6.18) of the order parameter M2{N), and by the exponential decay (6.25) 
of the stationary probabilities in all directions away from both endpoints. The mass, or 
inverse decay length, /i((/?) depends on the orientation: it decreases monotonically between 
(f = tt/A (decay along the diagonal j + k = N) and (p — (decay along both co-ordinate 
axes). 

This point brings some more insight on the flipping mechanism. Since the decay length 
is the largest along the axes, the walker, starting from (0,A^), will preferentially diffuse 
down the A;-axis, until it possibly reaches the vicinity of the origin, corresponding in the 
toy model to few (— ), few (+) and many holes in the system. Then, in a few steps, it can 
pass the line j = k and reach a region where the bias drives it to the j-axis. The typical 
duration of such a successful path, from (0, N) down to the j-axis, is presumably of order 
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N, indicating that, when the system flips, it does so on a very short time scale compared 
to the average duration of time it spends in a long lived state, i.e., near the two endpoints 
(0,Ar) or (TV, 0). 

This argument also provides an estimate of the average time between two flips (Tq), i.e., 
the average lifetime of the walker in one of the two regions near the endpoints. Indeed the 
probability that the walker reaches the vicinity of the origin is proportional to exp(— A?"//), 
with fj, = //(O), therefore (Tq) should be proportional to ex.p{NiJ,). 

We indeed find that the average time (Tq) between two flips in the toy model reads 

(To) « L Ar3/2 exp(iV//), (8.1) 



with 



H = //(O) = Inxi = In 



^(l-a- V(l-a)(l-2a)) 



(8.2) 



This result agrees with the estimate (1.4). It also confirms the above picture of the flipping 
mechanism. The mass /i can be interpreted as an activation energy, or barrier height per 
particle; the exponent 3/2 of the power-law prefactor is universal, i.e., it does not depend 
on the strength of the bias; the absolute prefactor L is only known in the continuum limit 
of a weak bias (a ^ 0), where L ^ 

We have thus demonstrated the occurrence of spontaneous symmetry breaking in the 
/? — > limit of the two-species exclusion model, in the sense that the system flips between 
two symmetry-related long lived states, by showing that the two properties mentioned 
above hold. 



Let us now impose boundary conditions which favor one species of particles over the 
other, breaking thus the symmetry explicitly. We consider an exclusion model where, 
instead of having boundary rates a, j3 for both types of particles, we have a, j3j^ for the 
positive particles and for the negative particles. We set (3± = H), where H is 

a symmetry-breaking field, such that \H\ < 1. Hence the positive (respectively, negative) 
particles are favored by the boundary conditions for H > (respectively, H < 0). If we 
now take — > 0, it can easily be checked that we find a toy model with the following 
hopping rates: a(l -|- H) west; a(l — H) south; 6(1 -|- H) north-west; 6(1 — H) south-east, 
again with a + b = 1/2. 

The bias associated with this model now reads 

v=«(:J)+(i-«)^(~i^)- (8.3) 

As long as the field is small enough, i.e., \H\ < Hq = a /{I — a), both components of the 
bias are still negative, and the mechanism which localizes the walker near the endpoints 
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is retained. The significance of this is as follows. Assume now > for definiteness. 
Even though the boundary conditions favor blocks of positive particles, the system may 
still spend lengths of time which diverge exponentially with the system size in a state 
dominated by negative particles. It can indeed be shown that the average durations of 
time spent in the (+) phase and in the (— ) phase scale respectively as (To)+ ~ exp(A'"//+) 
and (To)- ~ exp(A^/7,_), where fi.^. > /i_ are the masses along both axes, which are now 
different. This reinforces our proposed interpretation of the mass as the activation energy 
per particle, in the sense that Afi = fi^ — fi- is naturally interpreted as the difference 
in energy per particle between both pure phases. In the continuum limit of the model, 
defined as a ^ and H ^ simultaneously, we have fi± ~ {a±H)'^/ (4a), hence A^u ~ H. 

The analogy between H and a symmetry-breaking field in equilibrium statistical physics, 
such as a magnetic field in the Ising model, can be pursued. Consider the order parameter 
Mi{N) = {{k — i)/N), analogous to a magnetization. This quantity can be estimated by 
noticing that the statistical weights p± of both phases are proportional to the durations 
of time {Tq)± spent in them, hence Mi{N) pj^ — p_ k, tanh.{NH/2) within exponential 
accuracy, the latter formula being identical to the equation of state of the Ising model near 
coexistence. 

We finally notice that the unfavored phase is metastable for H < Hq and just unstable 
for H > Hq. In the latter case the system will always mostly consist of (+) particles. 

We close up with a short discussion on the behavior of the mass in the continuum limit 
of a weak bias. We consider the most general biased random walk in the quarter plane, 
where the velocity is still given by eq. (3.5), while the diffusion coefficients Dy and Dj^ 
are now different from each other a priori. It can be shown that the on-axis mass fi which 
enters the law (1.4) vanishes linearly with a, with a non universal prefactor depending on 
the diffusion coefficients, 

7 ^ -T- (8-4) 

2D^ (1 + ^1 + D||/D^) 

In the case of the toy model, with Dy ^ and ^ 1, we recover fi ^ a/4 in the 
continuum limit. Another case of interest is the north-south-east-west model, defined by 
moves to the four cardinal directions with rates b (north and west), a + b (south and east), 
(with a + 26 = 1/2), for which interesting properties can be derived by means of several 
approaches [15, 16]. In that case, we have Dy pa D_l pa 1/4, hence // ^ 2{\/2 — l)a in the 
continuum limit. 
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Appendix A 

In this Appendix we show that the toy model is the P ^ hmit of the two-species 
exclusion process, in a sense that will be made precise shortly. The argument implies that 
this will be true even if the time over which the models are compared, and the system size, 
become infinite as the limit is taken. 

We begin with a summary of the basic strategy. Let us write r]{t) for the configuration of 
the exclusion process at time t; r]{t) is determined by the initial configuration r]{0) together 
with some specification of the (random) transitions during the evolution. Similarly, X(t) = 
{j{T), k{T)) will denote the position of the walker in the toy model at time r, and is 
determined by X{0) and the random steps taken by the walker. Now to each configuration 
X — (j, k) of the toy model there corresponds a configuration ^{X) of the exclusion 

process: ^{X) = ( — 0---0 + h), with j (— ) particles and k (+) particles. We 

will show that T](t) agrees, with a probability which approaches 1 as /3 ^ 0, with the time- 
rescaled image $(X(2/3t)) of the toy process, assuming (a) that the processes coincide 
initially, i.e., 77(0) = $(X(0)), and (b) that the random transitions in the two models 
correspond. Achieving the latter is called coupling the models. 

It is convenient to realize the randomness in the exclusion process as follows. There are 
3(A'" — 1) + 4 distinct types of elementary transitions which occur in this process: three 
types of exchange (+ with — , + with 0, or with — ) for each nearest neighbor site pair, 
and entrance or exit of a particle at the left or right boundary. With each type of transition 
wc associate a "Poisson alarm clock," which is in fact just a random set of positive times, 
Poisson distributed with density equal to the rate — g, 1, a, or /3 — for the transition. The 
alarm clock "rings" at each time t* in this set, and the transition occurs if it is permitted in 
the configuration r]{t*) at that time. (Here and below we assume, for definiteness, that r]{t) 
and X{t) are left continuous, so that if a transition from configuration ( to configuration 
(' occurs at time t, then rj{t) = ( and r]{t+) = lims\tvi^) = CO We will refer to these 
random times generically as transition times and more specifically as entrance times, left 
exit times, etc.; as emphasized above they are in fact potential transition times. 

To couple the two models we use the random times from the exclusion process to deter- 
mine the evolution of the toy model; this coupling corresponds to the intuitive description 
in Section 2. Specifically, if t* is a left exit time, then r* = 2pt* is a time at which j can 
decrease in the toy model. Suppose that X{t*) — {j, k). There are several possibilities for 
X(r*+) = (/,/c'): 

Case I: j = 0. Nothing happens in the toy model: (/, k') — {j, k). 
Case II: j = 1. In this case, (/, k') = (0, N). 



34 



j 


Allowed transition at time 
t** in exclusion process 


k 


Final state 
(i(r*+),A;(r*+)) 


Case 


j = 


— 


N 


(0,iV) 


I 


i = i 







(0,iV) 


II 


i >2 


^ + 




{j-l,k + l) 


III.l 


^ - ••• 


k>0 




III.2.a 


k = 


(Ar,0) 


III.2.b 



Table A.l: Transitions in the toy model for which j may decrease. 
The initial state is {j,k) = {j (t*) , k{T*)) 



Case III: j >2. The transition in the toy model is determined by looking ahead of t* to 
the first transition time t** > t* for one of two types of transitions: either for the entrance 
of a + particle, in which case (j', A;') = (j — 1, A; + 1), or for a — — > — exchange on 
bond <1, 2>, in which case (/, k') = {j - 1, k), if k > 0, and (/, k') = {N, 0), if /c = 0. 

These possibilities are summarized in Table A.l. Right exit times correspond similarly to 
times of possible decrease of k in the toy model. Since the relative probabilities of cases 
III.l and III. 2 are a/{l + a) and 1/(1 + q;), this prescription reproduces the transition rules 
(2-4). Here it is important that transitions in the toy model are governed by the alarm 
clocks, not by the actual occurrence of transitions in the exclusion process, since the latter 
depend on the history of that process and hence their true probabilities are not known. 

Suppose now that we observe the two systems over a time interval [0, T] of the toy 
model— that is, for < r < T or < t < T/{2(3). Let B C [0, T] be the set of (bad) 
times T at which 77(r/(2/3)) 7^ $(X(r)), i.e., at which the systems do not coincide, i? is a 
random set since it depends on the transition times. Now we would not expect agreement 
at all times, since after the exit of a particle the exclusion system must rearrange itself 
before it will again agree with the toy configuration; we can thus at most hope that B 
have small Lebesgue measure Moreover, we cannot rule out completely the possibility 
that \B\ will be large, since if there are many exits in the exclusion process before the 
necessary rearrangements occur, the eventual configuration may not agree with that of 
the toy model. These considerations motivate the following notion of near- agreement: for 
e > we say that the two systems are e-close if the probability that \B\ > e is less than e. 

Now we can state the main result of this appendix. We consider a limiting process in 
which P approaches and admit the possibility that T = T{(3) and N = N{(3) increase 
during the limiting process; in particular, they may approach infinity, but not too rapidly. 
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Theorem: Suppose that r]{0) = ^{X{0)) and that T{(3) and N{(3) are monotonically 
increasing functions which satisfy lim^^o PNT = lim^_»o /?T^ = 0. Then for any e > the 
two systems are e-close for sufRciently small (3. 

Proof: For uniformity of notation it is convenient to assume that T(/3) — > oo as /? ^ 0. 
This entails no loss of generality, since we may if necessary replace T{j3) by T'{I3) > T{j3) 
with T'{(3) — > oo and with T'(/?) still satisfying the given hypotheses (this is the only place 
that monotonicity is used), and observe that if the two systems are e-close for T'{(3) then 
they are also e-close for T{(3). 

Let ti,t2,. . ■ be the exit times for the exclusion process, in increasing order, and let M 
be the number of these which occur between and T/{2(3). M is a random variable having 
Poisson distribution with mean and hence variance T, and Chebyshev's inequality implies 
that the probability that M < 2T is at least 1 — and hence goes to one as /? — > 0. Let 
us denote by Ei the event that M < 2T. 

Now let r denote the minimum of the rates 1, and a, and define 



Because 1 <^ a <^ /?~^, a time interval of length a is small on the r scale but large on 
the t scale. It is during such an interval that agreement between the two processes will be 
restored after a particle exit — that is, during an interval \ti,ti-\-a\. Restoration is simplest 
to analyze if no additional exit occurs during the interval; note that the probability that 
the next exit time tj+i satisfies tj+i < t j + a is 1 — e~^" < f3a. Thus the probability of 
the event that > + a for all i satisfying 1 < i < 2T is at least 1 — 2T/3a, which by 
hypothesis goes to one as /? — > 0. Let us call this event E2. When both Ei and E2 occur, 
all the relevant exit times differ by at least a. 

Suppose now that the two systems agree at the exit time tj, i.e., 'r]{ti) — $(X(Ti)), 
where = 2/3^^. If an exit in fact occurs in the exclusion process at time ti — for a left 
exit, in cases II and III of Table A.l — then agreement will not hold immediately thereafter. 
For each possible X{Ti) = {j,k), however, and knowing whether ti is a left or right exit 
time, we may identify a (j, k)- restoration event: one or more sequences of entrance and 
exchange times whose occurrence during the interval (ti, ti + a) will, if no exit time lies in 
this interval, guarantee the restoration of agreement by time ti + a. We will describe these 
restoration events below and prove that the probability that a (j, /c)-restoration event does 
not occur in a time interval (t, t + a) is, for small uniformly bounded by 2{ANI3 + 0^1"^). 
Thus the probability of the event that a X(Tj)-restoration event occurs in the interval 
(ti,tj + a) for all % satisfying 1 < z < 2T is at least 1 -4T(Ar/3 + /3V2)^ which by hypothesis 
goes to one as /5 — > 0. Let us call this event E's. 




(A.l) 
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We can now complete the proof. Take (3 so small that the probability that all of 
and £'3 occur is greater than 1 — e, and consider the evolution of the system when 
these three events do occur. Since the two systems agree at time zero they must agree 
until the first exit time ti. Since no exit, and an X(Ti)-restoration event, occur during the 
interval + a), agreement will be restored by time ti + a, and will persist until time 

^2- Continuing in this way, we see that agreement will hold for all t between and T/ (2(3) 
except for at most 2T time intervals of length a. Thus the total time of disagreement is 
at most 4f3aT on the r time scale. Finally, we take /3 so small that 4/3aT < e. 

It remains to consider the restoration events and estimate their probabilities. In the 
calculation we will assume that all entrance and exchange processes occur at rate since 
in so doing we can only underestimate the probability of restoration events. Again let 
ti — ''"i/(2/9) be an exit time — to be definite, a left exit time — and set X{Ti) — {j, k). We 
divide the discussion according to the cases of Table A.l. 

Case III: j >2. Agreement will certainly be restored if two sequences of transition times, 
corresponding to cases III.l and III. 2, both occur in the interval: (1) a left entrance time, 
then H — — > — h exchange times for bonds <1, 2>, . . . , <j — 1, j> in succession, then 

+ — > + exchange times for bonds <j, j + 1>, . . . , <N — k — 1, N — k>, and (2) 

— exchange times for bonds <1, 2>, . . . , <j — 1, j> in succession, followed if A; = (and 
hence j = N) hy a right entrance time. Each of these sequences requires at most N + 1 
transition times and hence has probability of nonoccurrence at most the probability that 
U < N +1, where U is a Poisson random variable with expectation (and hence variance) 
ra; by Chebyshev's inequality, 

rn AM 4- fi~^/'^ 

^ < ^ + ^ i ^ (.„ - (iv + 1))^ ^ = + 

Thus the event that one or the other sequence does not occur has probability at most 

Case II: J = 1 . Agreement will be restored when the system fills with N — k + particles. It 
is easy to enumerate the entrances and exchanges required, but these may occur in many 
different orders, so that a direct computation of the time needed for restoration appears 
difficult. We adopt an alternative approach by comparing the model to a well-understood 
system: a totally asymmetric one-species exclusion process on an infinite lattice, in which 
particles interchange with holes to their right at rate r. We suppose that sites in the 
one-species system are independently occupied or empty with probability 1/2 — this is a 
steady state — and couple the two systems by assuming that particle-hole exchanges on 
bonds <1, 2>, . . . , <A'" — 1, A'"> can occur in the one-species model only at -|-0 — > 0-|- 
exchange times for the same bonds in the two-species model, and that a -|-0 — > 0-|- 
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exchange on bond <0, 1> in the one-species model can occur only at a left entrance time 
for the two-species model. If we let J{u) denote the number of particles which actually 
enter on the left in the two-species model during the time interval {ti,ti + u), and J'{u) 
the number of exchanges which actually take place on bond <0, 1> in the one-species 
model during this same time interval, then it is easy to verify that J'{u) < J{u) as long 
as J{u) < N — k. From this it follows that the probability that J(a) < N — k, i.e., 
that restoration does not take place during the time interval {ti,ti + a), is at most the 
probability that J'{a) < N — k. But it is shown in [17] that J'{a) has expectation ra/4 
and variance V{a) satisfying lima^oo V{a)/ra — 0; thus again by Chebyshev's inequality, 

ProblJ(a) <N-t]< <^=^NP + (A.3) 

for j3 sufficiently small that V{a) < ra. 
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Appendix B 



This Appendix is devoted to an algorithm which allows for a fast recursive solution of 
the difference equations (4.10). 

• We start with the vj^k for the sake of simplicity. Setting j + k = n, and introducing the 
notation Vj^k = the recursion relation (4.10b) reads 

i^^('^) = Ff-'\ (B.l) 

with the boundary conditions Vq"^ = 1, vi"^ = 0. In eq. (B.l) D is the linear operator 
defined as 

Dvf=v^-^-b{v\%+v^-^,), (B.2) 

and the source term reads 

j-,(n— 1) f (n—l) , (n— 1)\ /-r, o\ 

F> ' = a IVj + ' j . (B.3) 

Eq. (B.l) is a linear equation for the Vj^\ with a source term F^""^) involving the 
vj^~^\ whence the possibility of a recursive solution. 
The solution to eq. (B.l) is the superposition 

) = {vnp + {vip , (B.4) 

of the solutions of the following two equations 

D{vn)'f^=0 with (^h)^"^ = 1, (^H)i"^ = 0, (B.5a) 
L» (^i)H = with (^i)^") = 0, {vit^ = 0. (B.5b) 

(with H for homogeneous, I for inhomogeneous) . 

The solution of the homogeneous equation (B.5a) goes as follows. Searching a solution of 
the form 2-^, we find z ~ exp(±a"), with the notation (6.12). Hence, imposing the boundary 
conditions, we obtain 

(vnp = (B.6) 

where 

= (B.7) 

smhcT 

is the Chebyshev polynomial of the second kind and of order n in the variable 

C = 2cosh(7= 1/6. (B.8) 
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These polynomials obey the recursion relation 



Un+l + Un-1 = CUn, (B.9) 



With UoiO - 0, UiiO = 1, SO that ^72(0 = C, Us{0 = C - 1, etc. 
The solution of the inhomogeneous equation (B.5b) reads 



n-l 

{vit^ = J2G^;;^Ft-'\ (B.IO) 

£=1 

where G^^) is the Green's function of D, defined by 

DG^;:}=5,,e with Gg=G3 = 0. (B.ll) 

The Green's function can be obtained explicitly as a superposition of plane waves exp(±j(T) 
in each of the sectors defined by the inequalities < j < £ and £ < j < n. We thus obtain 
the simple expression 

Gj,e = — Uini{j,£)Un-supij,i), (B.12) 

where mf{j,£) denotes the smaller of the two integers, i.e., £ if £ < j and j otherwise, 
whereas sup{j,£) denotes the larger, i.e., j ii £ < j and £ otherwise. 

The full recursive solution to eq. (B.l) thus reads 

(n) Un-j , COShcT-l^ / (^-1) , m 1 Q^ 

= -F7^ + rj t/inf(i,£)t^n-sup(i,£) \v\ + v\_^ > . (B.13) 

n n 

• Likewise one may write a recursion relation for the u^J^\ namely 



2 cosh ^ 



=1 

n-l 



1=1 



cosh cr — 1 



(B.14) 



One may deduce from the above analysis that the and are rational expressions 
in a, or equivalently in h. Their denominators are products of the Chebyshev polynomials 
up to Un- As a consequence, as n gets large, the poles of the vj"^ and become 
everywhere dense in the range — 1 < 1/(26) < 1. In terms of 6, the poles cover the whole 
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real axis except the interval [—1/2, 1/2]. In terms of a, the poles cover the whole real axis 
except the interval [0, 1] . 

The recursive solutions derived above provide a fast and efficient way of computing the 
Uj^k cind the Vj^k numerically. 

• The stationary probabilities p^^^ can be determined by means of a similar recursive 
scheme. First, eq. (6.1) is a closed-form equation for the p^^\ whose solution reads 



p(0)^^(0)Ci+^_ (B.15) 



obtain the following recursion relation for the p^'^^ 



Second, the structure of eq. (6.1) for u > 1 is very similar to that of eq. (B.l). We thus 

H cosher -1 / {u-l)\ /o i«\ 

Pj = fj-^ Ui^f{j,e)UN-u-sup{j,e) [Pi '+P}-i')- (B.16) 

,(0) ; 



=1 



Third, the first probability Pq is then determined by means of the normalization condition 
(4.5). 
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Figure captions 

Fig. 1. Characteristic curve D{x,y) = for a = 0.4 (see eq. (6.33)). The dashed hnes are 
the asymptotes. 

Fig. 2. Logarithmic plot of (To)a^ as a function of Na for the continuum hmit (full curve) 
and for the discrete random walk, with a = 0.1, 0.06, 0.03. 
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